Apparent violation of equipartition of energy in constrained dynamical systems 
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We propose a planar chain system, which is a simple mechanical system with a constraint. It 
is composed of masses connected by A'' — 1 light links. It can be considered as a model of a 
chain system, e.g., a polymer, in which each bond is replaced by a rigid link. The long time average 
of the kinetic energies of the masses in this model is numerically computed. It is found that the 
average kinetic energies of the masses are different and masses near the ends of the chain have large 
energies. We explain that this result is not in contradiction with the principle of equipartition. The 
apparent violation of equipartition is observed not only in the planar chain systems but also in other 
constrained systems. We derive an approximate expression for the average kinetic energy, which is 
in qualitative agreement with the numerical results. 
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PACS numbers: 

Information on energy distribution in many-body sys- 
tems is quite important for both theoretical and practical 
purposes. If a system is in thermal equilibrium, then, ac- 
cording to the principle of equipartition of energy, the 
average kinetic energy is equally distributed among all 
the degrees of freedom. Even when this principle holds, 
however, we found that nonuniform distribution of the 
average kinetic energy can occur. 

In this letter, we introduce a system called a "planar 
chain system" , which is a simplified model of a chain 
system e.g., a polymer. We show that the average ki- 
netic energy in this system is nonuniformly distributed 
even when it is in thermal equilibrium, but the principle 
of equipartition is not violated. We explain the reason 
for the nonuniform distribution of energy, which we refer 
to as the ^^apparent violation of equipartition of energy" . 
This property of apparent violation of equipartition of 
energy could provide a new insight into the behavior of 
chain systems. 

Let us introduce the planar chain system. The planar 
chain system is composed of N particles (masses) con- 
nected by — 1 links. The masses can rotate smoothly, 
as shown in FiglT] The links are massless and have 
fixed lengths. The system is defined by the following 
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Lagrangian L and constraints gi {i ~ 1,2, N — 1) 
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where N is the number of particles, TOj is the mass of I'th 
particle, Vi = {xi,yi) represents the position of the i'th 
particle, and £i is the length of the i'th link. U represents 
potential energy. We consider (i) a free chain with [/ = 
and (ii) external potential U = X]i=i ^(^)- 

If we define ipi as the angle between the j'th link and 
the —y direction (Fig[T]), we can rewrite the Lagrangian 
without the constraint. First we consider the following 
relations: Xi+i - Xi ^ li sin (p^ , y^+i - y., = -li cos . 
Using the total mass M and the center of mass [Xq, Yq) 



defined as M 



Xc 



i=l M 



Ei=i we obtain 

JV-l 

Xi = Xg + ^ aij sin ipj 
where a,-,- is defined as 
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i ifij , (3) 
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FIG. 1: A planar chain system 



By a straightforward calculation, we obtain the La- 
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grangian ([T]) in terms of (pi's and {Xg,Yg) as 
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-C/(Xg,1g, , (6) 

Ajki^p) = M|,i„(j, fc) Mmax(j,fe) cosiipjkyjh , (7) 

where cpj^ = ipj - ipk- 

We can consider this system as a simphfied prototype 
of various chain systems, e.g., proteins, polymers and 
spacecraft manipulators, under the assumption that the 
frequencies of bond-stretching vibrations are quite high. 

Now, we describe a method for numerical simulation. 
The Lagrangian that is expressed in terms of angles ([6]) 
is complicated and it is difhcult to numerically integrate 
the equation of motion, in particular for large N . Hence, 
we use the original form of the Lagrangian ([T|) and the 
constraint gi ([2]). Then, the equation of motion includes 
terms of the constraint, which is called a "Lagrange mul- 
tiplier" [l|. We determine Lagrange multipliers numer- 
ically at each integration step so that the constraint is 
satisfied 0. Methods of this type, e.g., "SHAKE" and 
"RATTLE" algorithms, are widely used for molecular 
simulation in chemistry [1, [3j 01 ■ In addition, some of 
the algorithms are known to be symplectic 0. Here, we 
use the forth-order symplectic integrator. In some cases, 
we verify the results by using an implicit Runge-Kutta 
method. 

If J7 = 0, the total angular momentum is conserved, 
hence, in this case, the energy distribution is different 
from the microcanonical distribution. In actual compu- 
tations, we place the system in a potential wall com- 
posed of arcs of radius 
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. Then, the system exhibits 

strongly chaotic motion similar to billiards 0] and does 
not have any conserved quantities other than the total 
energy, hence the microcanonical distribution is restored. 



between various parts of the system. Fig.[3]shows a power 
spectrum of xi{t) with the external potential mentioned 
above. It is a broad continuous spectrum, which is a 
manifestation of chaotic motion 
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FIG. 3: Power spectrum of xi{t) for < t < 32767. iV = 5. 
External potential in Fig[2]is used to obtain the spectrum. 

Using the method described above, we compute the 
long time average of kinetic energy. If the averaging time 
is sufficiently large, the long time average and thermal 
average can be assumed to be the same. 

The kinetic energy of i'th particle is defined as Ki{t) = 
^ (if + yf) , and its long time average is defined as 



K,{t)dt , t„ 



(8) 



K 



0.008 
I 0.007 
0.006 
0.005 
0.004 
0.003 
0.002 



♦♦♦♦♦♦♦♦♦♦♦♦ 



2 4 6 8 10 12 14 16 




FIG. 2: Planar chain in potential wall. N — 5, Af„aii = 4, 
a = 4Ne, R j = (i?,0),(-i?,0),(0,i?),(0,-i?), R ^ N£ + 
^a? - mP. £i = . . . = Im-i = L 

Although the planar chain system is a simple system, 
its dynamics is complex; further, energy exchanges occur 



FIG. 4: Long time average of kinetic energy Ki vs. i ((8|. 

= 16. rrii = 1.0 for all i and li = 1.0 for all i. The initial 
condition is as follows: Xi = (i — V) — N/2, yt — 0, p'f'' — 
for all i, p[^^ = -0.1, and p'f'^ = 0.1(i > 1). Here, ' and 
p\^^ represent the x and y components of the momentum of 
the i'th particle, respectively. The time step for integration is 
dt = 0.001. fmax ® is 10^. The relative error for total energy 
(square root of the time average of the squared displacement) 

is ^JaE'^/E^ = 6.1 X 10"". 

One might think that the values of all Ki^s in this sys- 
tem must be the same, by regarding Ki as the kinetic 
energy of the i-th. degree of freedom and applying the 
principle of equipartition of energy. However, this is not 
true. FigH] shows a plot of the average kinetic energy of 
each mass Ki ^ against i for = 16 planar chain sys- 
tem. It is clear that the Ki's are not equally distributed. 
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However, in the case of a planar chain system, equipar- 
tition of energy has a different meaning. From ([6]), we 
obtain the canonical momentum pi that is conjugate to 
t/jj as = ^ = J2k=i ^ikif)'fk , and we obtain the 

canonical kinetic energy A'j'^'' for the planar chain system 
as 
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FIG. 5; Convergence of kinetic energy Ki ([8]) in Fig|4] as a For example, eq. (|13p with = 3 and i = 1 we have 



function of t„ 



More importantly, we find that masses that are near the 
ends of the chain have large kinetic energies. We obtain 
this result for all the computed system sizes ( < 64 
)• Fig. [5]shows the convergence of /o"""' Ki{t')dt' as 
a function of ^max- The values shown in Fig [4] are well 
converged. 

However, this remarkable result is not in contradiction 
with the principle of cquipartition of energy. The prin- 
ciple of cquipartition of energy is stated as follows Q: 
Suppose we have a system defined by a Hamiltonian 



H{q,p) = K{q,p) + V{q), K{q,p) = 



1 



i,3 = l 



(9) 

where pi and qi are canonically conjugate to each other 
and Af is the total number of degrees of freedom. If it is in 
thermal equilibrium at temperature T, then the following 
relation holds: 



1 dK\ 1, ^ 

— n- ) — —KbT 



2^' dp J 



(10) 



(Summation over the index i is not taken in the left hand 
side.). The symbol (•••) represents thermal average at 
T, and is defined as 



(/(g,p)) ^ ly f{q,p)e-P"dT 



(11) 



for any function f{q,p). Here, dT is a volume element of 
phase space, Z is a partition function, and /3 = l/fc^T. 

Let us define the "canonical kinetic energy" K. 
the "linear kinetic energy" Ki as 



and 



1 dK 



(12) 



respectively. Here, cquipartition of energy means that 
the average values of K\ ' 's are equal at thermal equilib- 
rium. 

For systems such as gas models or lattice models. 



and ifi'^' K, : 



hence, the principle of 



cquipartition (jlOp simply means that {jniivf) = ^ksT 
which is a commonly used form of cquipartition of energy. 
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(14) 

It should be noted that Kl ' is defined by variables 
of every part of the system, whereas Ki is defined only 
by the i-th particle. In other words, canonical kinetic 
energy Kf^ is extended, whereas linear kinetic energy Ki 
is localized. Hence it is obvious that Kf^ ^ Ki. Since 

(c) 

Kl obeys cquipartition of energy, we can consider that 
Ki does not obey this principle. 
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FIG. 6: Long time average of canonical kinetic energy K!^ 
(|13p vs i. This plot is obtained by using the same data as 
that in FigH 

Fig. [5] shows the long time average of Kf for the 
same time series as that in Fig. 01 It is clearly shown 

that k'i'^^ 's take almost the same value for all i. That is, 
cquipartition of energy is realized. Non-equipartition of 
energy clearly shown in Fig[4]does not imply that princi- 
ple of cquipartition is violated. In other words, although 
the system obeys the principle of cquipartition of energy, 
the values of the linear kinetic energy are different at dif- 
ferent points in the system. We refer to the variation in 
(Ki) under thermal equilibrium as the ''^apparent viola- 
tion of equipartition of energy^'' . 

We derive the nonuniformity of the energy distribution 
by analytical calculation. By a straightforward calcula- 
tion, we obtain 



■mi 
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N-l 
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atjOik {cos (ipjk) (fjifik) ■ (15) 
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To evaluate the second term, we adopt the following 
approximations: 



(cos {(fijk) ipjipk) =0 for j ^ k 



(16) 
(17) 




(The matrix is included in exp{—f3H).) These ap- 
proximations indicate that each link in the chain rotates 
independently. Then, we obtain 

ksT M I 2 

(18) 

Details of the calculation will be shown elsewhere [9] . 

Eq.((T8]) shows that the average linear kinetic energy 
(Ki) varies from point to point. In other words, "appar- 
ent violation of equipartition of energy" occurs in this 
system. 

If all the masses are the same rrii — m, then we obtain 



ksT N 1 2 
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N-3 



E 



N-j 



This expression implies the following: 

(Xi) >{K2)>---< {Kn-i) < {Kn) 



(19) 



(20) 



It is clear that (Ki) is large at the ends of the chain 
and small at the center of the the chain: this result is 
in qualitative agreement with the result of the numerical 
computation shown in Figd) 

In this letter we have numerically shown that Ki = 
-^Wj does not obey the principle of equipartition of en- 
ergy for the planar chain system. Moreover Ki of parti- 
cles that are near both ends of the chain is large. The 
nonuniform distribution of the linear kinetic energy is 
qualitatively explained by analytical calculation. 

The apparent contradiction is due to the difference be- 



tween Ki 



^m^Vi M,iiK^ ^^i — 

is caused by the presence of the coordinate q {if for pla- 
nar chain systems) in the expression of the kinetic energy 
K, due to the existence of the constraint. Further, the 
same numerical time series show that the average values 
of K^'^'^ are equal, i.e., the system obeys the principle of 
equipartition. 

It is clear that there are other models in which the val- 
ues of the average kinetic energy are not equal. These 
models are systems with constraints, where the expres- 
sion of the kinetic energy includes coordinates. In fact, 
it has been found that the behavior of linear kinetic en- 
ergy in a multiple pendulum system is similar to that in 



? and K? 



-jPi 



dK 



This difference 



the planar chain system 
analysis elsewhere 11 1. 



[10| ; and we will report detailed 
In polymer science, the three- 
dimensional version of this model is known as a "freely 
jointed chain" 13, We expect that the behavior of 
the kinetic energy in the freely jointed chain will be sim- 
ilar to that in the planar chain system. 

We have shown that in the planar chain system, the 
energy at the ends of the chain is larger than that at 
the center. This result may be considered rather trivial, 
because it may seem that the end parts can be moved 
easily. However, even in thermal equilibrium, where all 
degrees of freedom have the same energy on average, the 
energy at the ends of the chain is large. 

This result would have important implications for the 
dynamics of chain systems such as molecules, proteins, 
polymers, and some artificial objects. For example, in 
polymer science, it is well known that atoms situated 
near the ends of the polymer chain have characteristic 
behavior called the "end effect" [l3]- Apparent violation 
of equipartition of energy we found in planar chain sys- 
tems can be closely related to the origin of the end effect 
of polymers. 
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